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ABSTRACT 

We present indications of rotation in the galaxy cluster A2107 by a method that 
searches for the maximum gradient in the velocity field in a flat disk-like model of 
a cluster. Galaxies from cumulative sub-samples containing more and more distant 
members from the cluster centre, are projected onto an axis passing through the 
centre and we apply a linear regression model on the projected distances x and the 
line-of-sight velocities V. The axis with the maximum linear correlation coefficient 
fmax — max [r(V, x)] defines the direction of the maximum velocity gradient, and 
consequently it presents the major axis of the apparently elliptical cluster. Because 
the effects of rotation are subtle, we put strong emphasis on the estimation of the 
uncertainties of the results by implementing different bootstrap techniques. We have 
found the rotational effects are more strongly expressed from distances 0.26-^0.54 Mpc 
from the cluster centre. The total virial mass of the cluster is (3.2 ± 0.6) x 10 14 M Q , 
while the virial mass, corrected for the rotation, is (2.8 ± 0.5) x 10 14 A / (q. 
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1 INTRODUCTION 

Velocity gradients suggestive of galaxy cluster rotation were 
found in several studies (e.g. Kalinkov 1968, Gregory 1975, 
! Gregory & Tifft 1976, Gregory & Thompson 1977, Materne 
& Hopp 1983, Materne 1984, Williams 1986, Oegerle & Hill 
1992, Sodre et al 1992, Biviano et al. 1996, Tovmassian 
2002). It is commonly accepted opinion that there is no 
evidence of rotation of the galaxy clusters, however, if the 
galaxy clusters rotate their dynamics and dynamical evo- 
lution would be different. The velocity dispersion profiles 
must be corrected for the rotation and the corresponding 
virial mass estimation will be different. Many theoretical 
constructions are based on the assumption of no rotation - 
e.g. the infall models and especially the theory of caustics 
(Regos & Geller 1989, Diaferio & Geller 1997, van Haarlem 
et al. 1993, Diaferio 1999, Rines et al. 2003). 

Here we make an another attempt to reveal rotation in 
the galaxy cluster A2107. We have chosen this particular 
cluster to illustrate our method, because it is well studied 
and indications of rotations were already found (Oegerle & 
Hill 1992). 

The structure of the paper is as follows: in Section 2 
we present our method, the different bootstrap techniques 
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that we use to estimate the uncertainties of the derived ro- 
tational parameters are presented in Section 3. In the next 
Section 4 we present the data for A2107 and in Section 5 
our results. We finish with discussion and conclusions (Sec- 
tion 6). In order to compare our results with the previ- 
ous studies we assume an Einstein-de Sitter cosmology and 
Ho = 100 km s _1 Mpc -1 . For this model at the cluster 
redshift z = 0.04100 << 1, 10' correspond to 358 kpc. 



2 METHOD 

We consider a flat, disk-like galaxy cluster with regions with 
nearly solid body rotation. The main idea is to find the axis 
of the maximum velocity gradient which, in the disk-like 
model, defines the major axis of the cluster. The minor axis 
is the axis of rotation. 

It is known the rotational effects are weak and their 
search would be successful if some cumulative technique is 
applied. Let us take sub-samples {i} of the cluster mem- 
ber galaxies, arranged by their projected distance from the 
cluster centre. Thus any sub-sample {i} contains the first 
i galaxies out to a given projected distance di. The sub- 
samples {i} are not independent. 

Let us introduce an axis passing through the cluster 
centre and rotating around it and we define the positional 
angle tp in the usual way (anticlockwise from N) as shown on 
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Figure 1. A schematic and idealistic representation of a disk- 
like cluster that indicates the axes and the angles used in the 
method. The symbols "+" and "-" are velocity deviations of the 
galaxies from the mean cluster velocity. Note the positional angle 
of rotation (pf, spans the whole range from to 360 degrees. 



Fig. 1. We assume the axis has direction in sense SN to WE, 
NS, EW. In this case the positional angle is 0° < ip < 360°. 

We project the galaxies from sub-sample {i} onto an 
axis with positional angle ip and examine the linear regres- 
sion 

V = a + p{x-(x)), (1) 

where V is the observed line-of-sight velocity, x are the pro- 
jected cluster-centric distances on the chosen axis, (x) is the 
mean of Xi for the sub-sample {i}, and in fact a = (V). 

Our aim is to find the maximum gradient in the veloc- 
ity field. In linear approximation, this corresponds to the 
axis with iprmax for which the linear correlation coefficient 
between V and x is maximal, r max - Consequently, the stan- 
dard deviation (st.dev.) of the regression for the sub-sample 

s=(^p(Vk-{Vi)) 2 /(i-2) ] J (2) 

is minimal and we denote it by s rmax which corresponds to 
Tmax for the same ip; 14 is the line-of-sight velocity estimate, 
derived from the regression. 

If the velocity field of a flat cluster is influenced by ro- 
tation then the maximum velocity gradient will be along the 
major axis a (or equally, along the axis of maximum elonga- 
tion). The positional angle of the maximum velocity gradient 

iS ifirmax — ifia — (fismin , 

where the indexes rmax and smin 
denote the maximum correlation coefficient and the min- 
imum standard deviation of the regression estimate. Cor- 
respondingly, the axis of rotation is at <p r o = fib = fismax, 
where rO denotes zero correlation coefficient. Note that there 
are two positional angles where r is zero. We define the po- 
sition of the minor axis with ip r o > ip r max- 

Searching for the major axis or maximum elongation 
axis with a limited number of cluster members is prone to 
large uncertainties and few interlopers could bias the esti- 
mated ellipticity and shape. Much better approach would be 
by using X-ray observations because the hot gas maps bet- 
ter the cluster potential and subsequently the cluster shape. 



Unfortunately the positional angles of the maximum elon- 
gation in the cluster cores are known to vary depending on 
the X-ray isophote level. 

The virial mass of a cluster is defined through the 
velocity dispersion cry. In fact, 



av = s r0 ((i~2)/(i-l)) 1/2 



(3) 



If the cluster rotates the velocity dispersion has to be cor- 
rected (ay), according to the formula above, but with s rm ax 
instead of s r o- Then M. — M c = AM is a fictitious mass, 
due to rotation. 

Supposing that some noticeable effects of rotation re- 
ally exist, then, with our approach, we expect the following 
indications: 

(i) The correlation coefficient r max will be significant for 
at least few consecutive sub-samples {i}. 

(ii) The positional angles (firmax for these consecutive 
sub-samples will not be randomly distributed along the en- 
tire range of ip, but in a relatively narrow interval. 

(iii) The variation of r and s for 0° ^ ip < 360° would be 
close to sine waves and, for sub-samples {i} with significant 
r ma x, the phase shifts of both curves should be according to 
(ii), in a narrow range. The expressions for the sine waves 
of r and s are 



r(ip) = r max sin(ip - ip rma x) 
and 



(4) 



s(tp) = (sro + s rmax )/2 + [(sro - Wi)/2] sin(v3 - ip r o) (5) 

(iv) The virial mass for some sub-samples {i} will be sig- 
nificantly different as computed for s r o and s m „. 



3 UNCERTAINTIES, BOOTSTRAP 

TECHNIQUES AND RANDOMISATION. 

The question of uncertainties is a crucial one because the 
effects of rotation are weak. Where it is possible we rely on 
bootstrap uncertainties. 

We apply two bootstrap techniques. The first one is the 
standard or the classical bootstrap (B) - (Efron 1982, Efron 
& Tibshirani 1986). Shortly, for each sub-sample we find the 
positional angle ip rm ax and then resample the velocities V 
and the projected distances x on this axis in order to derive 
the bootstrap uncertainties of a, j3, r, s. This procedure 
however has the disadvantage that it is impossible to find 
the uncertainties of ip rma x and ip r o- 

That is why we apply simultaneously a modified boot- 
strap (MB) techniques: the positions and velocities of sub- 
sample {i} are subject to bootstrap resampling after which 
the search for r ma x is carried out and the corresponding 
firmax, Q, 13, r, s,... are derived. Note that each bootstrap 
sample leads to completely different r max , (firmax, a, [3,... 

In all cases, we run 10 4 bootstrap generations and cal- 
culate the corresponding quantities for ip[0° — 360°) with a 
decrement of 0?1. We use the resulting distribution of the 
quantity in question to derive the confidence interval. 

Let us denote with wb whatever of the quantities 
we discuss. The first estimator of the uncertainty we use 
is the bias-corrected (be) 68% confidence interval w £ 
[W6c(0.16), w bc (0.84)] (Efron & Tibshirani 1986) with 
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w bc (t) = CT 1 {$ + 2$- 1 (G((w o6s »)] } , (6) 

where G is the cumulative distribution function (CDF) of 
the Wb, {w bs} is the observed mean value, and $ is the CDF 
of the normal distribution (thus ^ 1 (t) = — 1, +1 for t = 
0.16 and t = 0.84 respectively). This uncertainty is Aiobt- 
Efron & Tibshirani (1986) uncertainty is effectively used by 
Shepherd et al. (1997) and Kalinkov, Valtchanov & Kuneva 
(1998a) for the space correlation functions of galaxies and 
clusters of galaxies. 

The second estimator of the uncertainty we implement 
is the bootstrap standard deviation (Ling, Barrow & Frcnk 
1986) 



Aw i 



1/2 



(7) 



where (wb) is the mean bootstrap value. This st.dev. is very 
effectively used by Kalinkov & Kuneva (1986). 

We have found that both uncertainties Awb and Auibt 
are consistent with each other, but we prefer Aw et, since 
the estimator is valid even if the distribution of w is not 
Gaussian. 

The uncertainty on the correlation coefficient can also 
be computed using the classical Fisher z-transformation (e.g. 
Press et al. 1992, Section 14.5). The uncertainties derived by 
this method are also consistent with the other two. 

There is one problem - the bootstrap technique does not 
work when we determine the uncertainties of the virial mass 
estimates. In the resampling, some objects are duplicated 
and consequently the distance between them is zero and the 
cluster potential tends to infinity. In this case only we apply 
the jackknife estimate (see e.g. Efron 1982). 

In order to assess the significance of the derived quanti- 
ties we use randomised sub-samples for which the azimuthal 
angles (or ip) of the galaxies centres are made randomly dis- 
tributed in the interval [0° , 360°). This new sample is subject 
to the same procedure of parameter and confidence interval 
estimates. Here the uncertainties are only Auib. It is worth 
to note that with this procedure the observed radial density 
and velocity distributions are not altered. 



4 DATA 

Coordinates and heliocentric velocities of galaxies in the di- 
rection of A2107 are taken predominantly from the most 
complete optical study of Oegerle & Hill (1992) - further on 
OH. Some corrections and additions are made according to 
NED (new redshifts and coordinates) and therefore new ob- 
servational errors are defined. We use also data from other 
sources, e.g. Zabludoff et al.(1993). We reckon that galaxy 
Nr. 252 is the same as Nr. 273 (LEDA 94259) from the list 
of Oegerle & Hill (1992). 

We adopt the X-ray centre of A2107 according to Ebel- 
ing et al. (1996) RA = 15 h 39 m 38!0 and Dec = 21°47'20" for 
J2000. Our results however do not substantially differ nei- 
ther if we adopt the original cluster centre of Abell (1958) 
and Abell, Corwin & Olowin (1989) nor if we take the cD 
galaxy (UGC 9958) as centre. 

The member galaxies were selected using ROSTAT 



(Beers, Flynn & Gerhardt 1990) as well as the cone dia- 
grams and the velocity distribution. There are 70 galaxies 
out to 37.0 arcmin from the X-ray centre within (V) ± 2<r„. 
Zwicky galaxy 136-22 is the farthest member. For the non 
relativistic case (i.e. in the observed rest frame), the mean 
velocity is (V) = 122911^5 km s _1 , the standard deviation 
627^4® km s _1 , the bi-weighted estimate of location and 
scale Vbw = 12318 km s _1 and Sb = 639 km s _1 , the median 
Vm = 12380 km s _1 and the normalised absolute deviation 
s m — 623 km s _1 , skewness -0.28 and kurtosis -0.40. 

All our results cited here are for this sample with size 
n = 70. 

There are two galaxies at 1.2 arcmin from the cluster 
centre (Nrs. 242 and 289 from the list of OH) with velocities 
14069 and 14028 km s _1 . These galaxies could be taken 
as cluster members if we stand on the infall theory, since 
the velocity dispersion at the cluster centre should be the 
largest. Including these two galaxies we obtain for n = 72: 
(V) = 12339lgf, km s _1 , s = 683±g, V bw = 12338, s b = 685, 
Vm = 12386, s m = 635 km s _1 , skewness 0.05, kurtosis -0.03. 
All further conclusions are also referred to the sample with 
size n = 72. 



5 RESULTS 

In Fig. 2 we show r max results for the sub-samples 6 ^ 
{i} $J 70. For sample {70} the correlation coefficients are 
r ma x,ro = 0.3581 am , rma X ,rand = 0.154 ± 0.077. Assuming 
that correlation coefficients r ma x,70 and r m ax,rand are almost 
normally distributed and are independent, the significance 
of their difference {r m ax) — (r ra nd) can be estimated from 
Student's distribution. Consequently, the probability of the 
null hypothesis, that the correlation coefficients are drawn 
from one and the same population, is V < 1.0 ■ 10~ 13 . 

The highest r ma x,47 = 0.637lg'joo (where st.dev. are 
according to the MB method) occurs for sub-sample 47, at 
di = 0.54 Mpc. For B the st.dev. are (+0.069,-0.084). For 
the randomised sub-sample 47 we have r m ax,rand = 0.209 ± 
0.104. In this case V < 1.0 • 10~ 17 . 

A crucial test for rotation are the diagrams {ip r max, di) 
and (ip r o,di). If the points are randomly distributed in 
the interval [0° , 360° ) the hypothesis for rotation must be 
rejected. The results for both (p r max and (p r o are shown 
on Fig. 3, together with the corresponding results for az- 
imuthally randomised sub-samples (denoted as stars). Note 
that ip r o,rand = <Prmax,rand- The st.dev. for the randomised 
sub-samples are « [(360°) 2 /12] 1 / 2 which corresponds to the 
st.dev. of uniformly distributed positional angles. The dif- 
ferences between ip r Q and ip r max are « 90°. 

Next, for sub-samples with 0.54 < di Mpc < 0.76, for 
which r ma x significantly differs from the azimuthally ran- 
domised correlation coefficients, (tp r max) = 173° ±4°. We as- 
sume that this is the positional angle of the velocity gradient. 
The positional angle of the minor axis is {ip r o) = 271° ± 5°. 
The last value gives the positional angle of the rotational 
axis if our model is valid. The difference between both mean 
values is 98°. 

In Fig. 4 we show the st.dev. of the regressions s rm ax 
and Sro for positional angles corresponding to r ma x and rO 
- there is significant distinction for those sub-samples where 
the maximum correlation coefficient is significant. 
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Figure 2. Observed r ma x (dots) and azimuthally randomised 
( r max,B) (diamonds) for subsamples {i} = 6 4- 70 as a function 
of the distance from the cluster centre di. The st.dev. for r max 
are defined according to MB while for (r max g) are according to 
B. These types of errors are presented on all figures except Figs. 
6-7. 
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Figure 3. The positional angles ip rm ax (dots) and ip r o (dia- 
monds). The stars show the results for azimuthally randomised 
sub-samples with the errors (dotted lines). 
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Figure 4. St.dev. of the regressions: s r o (dots) and s r 
monds). 
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Figure 5. Variation of r (solid line) and s (dashed) as a function 
of the positional angle ip for sub-sample {i} = 47. 
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Figure 6. The virial masses: A4 C (dots) and A4 (diamonds). The 
error bars refer to the jackknife estimator. Fits with third order 
polynomials are also shown. 



The variations of r and s as a function of ip for sub- 
sample {i} = 47 are given on Fig. 5. The curves are very 
close to sine waves. For the maximum correlation coeffi- 
cient determined from azimuthally randomization of any 
sub-sample {i}, the s curve is also sine- like but with much 
smaller amplitude and for various bootstrap resamples their 
phase angles are randomly distributed on [0°,360°). 

For the sub-sample {47} we have <p rm ax = 166° and 
(fro = 260°. And the virial mass estimate M = (2.18 ± 
0.36) • 1O 14 X and M c = (1.29 ± 0.21) ■ 10 14 A^ Q . 

We have computed the cumulative virial mass of A2107 
for sub-samples {20} to {70}, following the prescriptions 
of Heisler, Tremaine & Bahcall (1985). We have used the 
velocity dispersion ay to derive the usual virial mass M. as 
well as the corrected virial mass A4 C . The results are given 
on Fig. 6. The error bars were calculated implementing the 
jackknife method. There is a slight difference between M 
and M c . The distinction begins at {i} — 20, at d,2o = 0.22 
Mpc. In order to compare both estimates we fit the masses 
with third order polynomials, as shown on Fig. 6. On Fig. 7 
we show M — M c together with the difference between both 
polynomials. If our finding that the cluster A2107 rotates 
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Figure 7. The difference between M and A4 C . 
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Figure 8. Parameters /3 rma x (dots) and (Prmax, mb) (dia- 
monds) . 



is true, then Fig. 6 would reflect a fictitious mass, due to 
rotation. Our estimates for the total masses are M. = (3.21± 
0.62) ■ 10 14 Af Q and M c = (2.80 ± 0.54) ■ 10 14 Af Q . 

The parameter (3 of the examined regression, which rep- 
resents the maximum velocity gradient expressed in km s _1 
Mpc -1 is shown on Fig. 8. Actually it is f3 rma x and also 
(Prmax,MB) as functions of di. 

Our method gives an opportunity to estimate an up- 
per limit for the rotational period. In Fig. 9 we show the 
upper limits for P rm ax compared with the results from az- 
imuthally randomised samples. The estimate is simplistic - 
we have not included any correction for the inclination of 
the galaxy cluster. Any inclination of our flat model will 
inevitably decrease the period. 

The analysis up to now was based on cumulative dis- 
tributions. It is possible however to obtain an additional 
information by considering the differential case. Let us de- 
note the differential sub-samples by {p,p + q — 1}, where 
p = 1,2, ...,70 and q = 3, 4, . . . , 70. Repeating our pro- 
cedure for the differential sub-samples, which are in fact 
galaxies in different annuli, we find that r max is very noisy 
for small p and q and at its maximum rmax = 0.8141^102 
for {25, 47} (i.e. p = 25 and q — 23 or for distances 
0.26 -T- 0.54 Mpc from the cluster centre). For the same sub- 
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Figure 9. Upper limits of the period P rm ax (dots) and 
(Prmax, mb) (diamonds). 



sample the azimuthally randomised bootstrap resampling 
gives r ma x = 0.271 ± 0.133. The corresponding probability 
is V < 1.0 ■ 10~ n and it is larger than that for sub-sample 
{47} since there are fewer degrees of freedom. 

The standard deviations of the estimates for the r max 
and ro are s m „ = 44614J 8 km s _1 and s r o = 767t|g 
km s _1 correspondingly. For the differential sub-sample 
^Prmax — 167° and tp r o = 261°. For these 23 galaxies the de- 
rived period is P = (2.431°;^ )-10 9 years. The corresponding 
differential virial masses are M = (3.51 ± 0.60) • 10 14 .Mq 
and M c = (1.19 ± 0.20) • 10 14 X Q . 

The results from the differential treatment are in agree- 
ment and support the cumulative ones. 



6 DISCUSSION AND CONCLUSIONS 

A2107 is a nearby cluster of richness class 1, classified 
as BMI (Abell 1958, Abell, Corwin & Olowin 1989). The 
Abell's count of galaxies is Na = 51, the cluster radius is 45 
arcmin (22.4 arcmin from Struble & Rood 1987, Ho = 50) 
and the RS classification type is cD (Struble & Rood 1987). 
The cluster is isolated and is also identified (Kalinkov, 
Valtchanov & Kuneva 1998b) as Zwicky cluster ZC 7573 
(Zwicky & Herzog, 1963) as distance group Near, type mc, 
population of galaxies 293 and equivalent radius r e z = 69.4 
arcmin. At 26.5' and positional angle of 204° is located an- 
other Zwicky cluster ZC 7578 (distance group VD, popula- 
tion 110, type mc, r e z = 7.8 arcmin). 

There is no detailed description of A2107. Girardi et al. 
(1997) comment on A2107 - "Remarkably regular cluster". 

Presumably rotation of clusters of galaxies is difficult 
to be established with certainty. Materne & Hopp (1983) 
have shown that it is extremely hard (if generally possi- 
ble) to distinguish the case of a single cluster in rotation 
from the case of two overlapping clusters, which are merg- 
ing or departing from each other. That is why OH, investi- 
gating the same cluster A2107, have examined both possi- 
bilities. They employed the <5-test of Dressier & Shectman 
(1988) and they have found that there is significant sub- 
structuring in A2107. Supposing that Ym=i ^ ~ n > wnen 
there is no spatial-velocity correlation, OH performed 1000 
Monte Carlo (MC) simulations and found (J"^ 8i)mc = HI 
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(see also Oegerle & Hill 2001). Because their sample con- 
tains n = 68 cluster members, they concluded that A2107 
contains substructures. All ^<5;,mc are smaller than 111 
and hence the probability for existence of substructures in 
A2107 is > 0.999. We have carried out exactly the same ex- 
periment but with 10 4 Monte Carlo simulations and found 
that ~^2 t 8i,mc are almost normally distributed with mean 
(J2i S ')mc = 70.00 ± 10.82 against X)<5i, oi)s = 102.3. Our 
result support the conclusion of OH that A2107 indeed con- 
tains substructures with probability > 0.994. 

Concerning the velocity distribution, relying on the 
large peculiar velocity of the central cD galaxy (UGC 9958, 
V pec = 270 km s" 1 ; see also Oegerle & Hill, 2001), OH de- 
composed the velocity histogram into two sub-clusters. They 
have chosen one decomposition amongst infinite number of 
possibilities. Indeed, there are no reasonable arguments that 
constrain us to assume two sub-cluster centres or that the 
velocity histogram is superposition of two (or may be more?) 
distributions. 

We have verified the hypothesis for sub-clustering - one 
against two centres. We have used another powerful test 
especially for two sub-structures (Lee 1979), introduced in 
the cluster sub-structure studies by Fitchett (1988) and suc- 
cessfully applied by Fitchett & Webster (1987), Rhee, van 
Haarlem & Katgert (1991b), Pinkney et al. (1996). Our re- 
sult, on the base of 10 000 simulations of azimutally ran- 
domised samples, leads to Lee function L az = 1.42 ± 0.21 
while Lobs = 2.111q'4^. The alternative hypothesis, that the 
cluster has two sub-structures, is accepted at significance of 
at least 95%. The corresponding positional angle of the two 
sub-cluster centres is 145°. 

But the subclustering does not rule out the hypothesis 
of rotation. It is plausible that the cause of sub-clustering in 
the Dressler-Shectman diagram (Fig. 5 of OH) is just due to 
the rotation, since the most prominent dumpiness is along 
the maximum gradient in the velocity field of A2107, indi- 
cated by our ip r max- Indeed, from galaxy studies we know 
that the velocity field of disk-like rotating galaxies without 
any massive halo has two extremes, V m ax and V m in- This 
would cause apparent sub-clustering. But if there are two 
overlapping clusters, then it will be unlikely they will gener- 
ate a mimicry of rotation - first of all, the velocity histogram 
will have two peaks or quite broad velocity distribution with 
rather unrealistically high velocity dispersion. None of this 
is observed in A2107. Secondly, we have investigated the 
behaviour of correlation coefficients, positional angles etc. 
From our point of view, if there are two sub-clusters that 
are located at positional angle close to <p rm ax, this would 
be regarded as evidence for rotation. Nevertheless some au- 
thors disregard this possibility as the cause of the velocity 
gradient (den Hartog & Katgert 1996). 

Our strongest evidence for rotation is the consistent 
measure of the positional angle of the velocity gradient 
<prmax for consecutive cumulative sub-samples. One way of 
validating the inferred positional angle is to compare it with 
the positional angle of the elongation of the cluster tp e i . For 
our flat disk- like model y> ei = ip r max = fa- Debating the 
Binggeli's effect (Binggeli 1982) many authors have mea- 
sured ip e i , as well as ellipticities, of various samples of Abell 
clusters (e.g. Struble & Peebles 1985, Struble & Rood 1987, 
Lambas et al. 1990). While Struble & Peebles (1985) have 
found <p e i = 137° for A2107, Rhee & Katgert (1987) quote 



160° and 93° for the 50 and 100 brightest galaxies respec- 
tively, within a radius of 1 Mpc. According to Binggeli (1982) 
(fei = 62° . In a detailed investigation Rhee, van Haarlem & 
Katgert (1991a, 1992) apply three different methods and de- 
termine for the same cluster, using about 300 galaxies out to 
0.75 Mpc, that ip e i = 75, 68 and 90° with ellipticities 0.13, 
0.00 and 1.00 respectively. But according to their Fig. 5 
(Rhee, van Haarlem & Katgert 1991a) <p e i « 143° for the 
first 20 brightest galaxies while for the first 70 is about 3°. 

Apparently the estimates of <p e i for A2107 are not very 
reliable. Perhaps the reason is objective - it seems that (p e i 
depends on the brightness of the galaxies and the distance to 
the cluster centre. There are examples of drastic variations 
of ip e i in some clusters in the most elaborate paper on this 
subject (Burgett et al. 2004). 

As stated above there are two possibilities to explain 
the velocity gradient in a galaxy cluster - two-body model 
or rotation. An anonymous referee however turned our at- 
tention on a third possibility: according to simulations of 
galactic tides in dwarf spheroidal galaxies Piatek & Pryor 
(1995) have shown that tides produce large ordered motions 
which induce apparent rotation. Biviano et al. (1996) as- 
sume that the velocity gradient in the central part of Coma 
cluster maybe due to tidal effects caused by falling groups 
of galaxies and not by rotation. But in the case of A2107 it 
is not quite the same, because we did not find any groups of 
galaxies in this regular cluster. 

The parameters j3 and <p rma x derived with our method 
may be compared with those of other authors. The positional 
angle of the rotational axis of OH, namely 70°, is close to 
our determination 271° - 180° = 91°. From dV/dX, dV/dY 
for A2107 in Table 3 of den Hartog & Katgert (1996) we 
infer f rma x ~ 330°. The last value is close to our estimate 
of (frmax = 173° with 180° difference depending on the con- 
vention used. 

The parameter f3 after den Hartog & Katgert (1996) 
is 842 km s" 1 Mpc" 1 (for H = 50 km s" 1 Mpc" 1 ). Our 
value (3 for the entire cluster A2107 is 718 while for the 
subsamples {47} and {25, 47} the values are 2232 and 2527 
km s" 1 Mpc" 1 correspondingly. OH give « 1800 km s" 1 
Mpc -1 for the whole cluster, which is close to our value for 
the sub-samples but not for sample {70}. 

Our estimate of the rotational period is P = (2.4^2) ' 
10 9 years, while the azimuthally randomised bootstrap gen- 
erations give P = (1.3 ± 2.3) • 10 10 years. We have to note 
that our method of searching the maximum correlation in 
random samples leads to artificially decreasing the period of 
rotation. The period quoted by OH is 7 • 10 9 years. These 
periods are not corrected for the inclination of the cluster. 
Our estimate of P could be compared with the correspond- 
ing value of Gregory & Tifft (1976) for the outer regions in 
Coma cluster (A1656), which is about 2- 10 11 years. We pre- 
sume the real period is < 2 • 10 9 years and therefore A2107 
has made already few rotations during the Hubble time. 

There are several virial mass estimates for A2107. All 
estimates are not substantially different because they are 
based on one and the same radial velocity data. According 
to OH the total virial mass of A2107 is M ~ 4 ■ 10 14 X Q 
and in the frame of their two-body model, the masses of 
both sub-structures are ~ 1.3 : 1. Girardi et al. (1998) give 
a virial mass M = (3.02l 9 'gg) • lO 14 ^© and a corrected 
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M c = (2.62±g;^) • 1O 14 X (according to their Eqn. 8). 
Girardi et al. (2002) use the corrected mass only. But the 
correction of Girardi et al. (1998, 2002) has quite a different 
nature - it is due not to rotation but to a modification of the 
virial theorem in order to compensate the incompleteness of 
the sample out to the virial radius. 

Our estimated total mass of A2107 is Ml = (3.21 ± 
0.62) ■ 1O 14 .M0 (without any correction). Most importantly, 
reducing the effect of rotation we infer virial mass M c — 
(2.80 ± 0.54) • 10 14 X Q . The simplest way to get this virial 
mass is to estimate the velocity dispersion from the linear 
regression at positional angle <p r Q, 
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